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Abstract 


Permeability  and  bulk  densities  of  duff  and  humus  in  stands  of 
ponderosa  pine,  lodgepole  pine,  and  Douglas-fir  were  measured  to 
provide  a  data  base  for  fuel  moisture  and  soil  infiltration  models. 
Diffusion  of  heat  and  water  vapor  through  a  ponderosa  pine  duff 
layer  was  measured  to  calibrate  a  diffusion  model.  These  results 
provide  necessary  data  and  relationships  required  by  complex  forest 
soil  water  balance  and  energy  models. 
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Introduction 


Heat  and  water  movement  in  forest  soils  are 
integral  parts  of  many  forestry-oriented  models. 
Knowledge  of  infiltration  is  a  basic  requirement 
of  hydrologic  and  ecological  prediction  models,  al- 
though soil  moisture  also  is  important  in  engi- 
neering through  analysis  of  soil  stability.  Influ- 
ence of  soil  moisture,  primarily  the  upper  hori- 
zons (L,  F  and  H),  on  fire  danger  and  fire  manage- 
ment require  detailed  and  accurate  prediction. 
While  fire  receives  the  major  emphasis  in  this 
paper,  applications  to  other  forest  management 
activities  also  are  included. 

Prediction  of  heat  and  mass  flux  in  forest  soils 
with  models  based  on  physical  processes  almost 
always  requires  soil  structural  data  which  can  not 
be  measured  readily.  The  most  practical  pro- 
cedure to  develop  this  necessary  data  base  would 
be  to  develop  similitude  relationships  between 
the  easily  measured  characteristics  of  depth  and 
bulk  density  and  those  properties  of  timelag,  per- 
meability, and  thermal  conductivity  required  in 
the  fluid  and  thermal  flux  models. 

Many  bulk  density  measurements  have  been 
made  for  litter,  duff,  and  humus  under  most  of 
the  major  conifer  species  and  some  of  the  de- 
ciduous species  (Van  Wagner  1970;  Stocks  1970; 
Brown  1966,  1970a;  Anderson  1969;  Davis  et  al. 
1968;  Ffolliott  et  al.  1968).  Laboratory  measure- 
ments of  particle  characteristics  such  as  surface 
area  to  volume  ratios,  moisture  and  thermal  time- 
lags  and  particle  densities  (Brown  1970a,  1970b; 
Mutch  and  Gastineau  1970;  Fosberg  et  al.  1970; 
USD  A  Forest  Service  1966;  Anderson  et  al.  1976; 
Blackmarr  1972;  Van  Wagner  1969)  have  been 
made.  In  addition,  several  studies  have  been  car- 
ried out  on  flux  properties  of  forest  soils  or  on 
sorbing  porous  media  which  describe  heat  and 
moisture  flux  of  specific  litter,  duff,  or  humus 
layers  (Nelson  1969;  Anderson  et  al.  1976;  Van 
Wagner  1970;  Johnson  1968;  Jarvis  and  Tucker 
1968;  McLeoud  1976;  Murano  and  Lawson  1970). 


A  few  studies  have  included  measurements  of 
parameters  such  as  the  void  volume  per  amount 
of  surface  area  (Anderson  1969;  Brown  1970a) 
which  can  be  directly  related  to  the  mass  and 
energy  flux  parameters.  These  data  from  the  liter- 
ature along  with  data  collected  in  this  study  were 
combined  to  develop  methods  to  calculate  time- 
lag  and  infiltration  parameters  required  by  heat 
and  mass  flux  models. 


Theoretical  Basis  for 
Data  Relationships 

Three  types  of  relationships  can  be  defined  to 
interrelate  variables  or  groups  of  variables.  The 
first  type,  called  state  relationships,  is  based  on 
the  fundamental  definition  of  the  individual 
parameters.  A  second  type,  defined  as  dynamic 
relationships,  relate  flux  parameters  to  the  indi- 
vidual parameters.  The  third  type  introduces 
time-dependent  processes  and  relates  the  state 
and  dynamic  relationships  to  time-dependent 
characteristics  such  as  timelag.  This  type  will  be 
called  similitude  relationships. 


State  Relationships 

Five  basic  particle  and  horizon  characteristics 
defining  the  soil  structure  are  required  to  de- 
scribe the  state  relationships  and  to  calculate 
combinational  state  characteristics.  These  struc- 
tural characteristics  are  particle  surface  area  to 
volume  ratio  (a),  particle  density  (@p),  porosity  (</>), 
bulk  density  of  the  horizon  (eB),  and  thickness  of 
a  horizon  (H).  The  first  two  characteristics  define 
the  particle  properties  and  are  readily  measured 
in  the  laboratory.  In  forest  soils,  these  properties 
in  the  litter,  duff,  and  humus  layers  are  fairly  con- 
stant and  can  be  related  to  the  plants  on  the  site. 
The  last  three  characteristics  define  the  way 
these  particles  are  assembled  in  a  soil  horizon  and 
are  measured  in  the  field. 
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Procedures  and  techniques  for  measuring  these 
properties  are  described  elsewhere  (Brown  1970a, 
1970b;  Siau  1971).  A  complete  set  of  definitions 
and  elementary  manipulations  are  given  in  the 
Appendix. 

State  relationships  require  detailed  laboratory 
measurements  of  particle  density  [qp)  and  surface 
area  to  volume  ratio  (a).  Because  these  character- 
istics are  well  known  and  are  reasonably  con- 
stant, only  field  measurements  of  bulk  density 
and  horizon  thickness  are  required.  Since  bulk 
density  measurements  for  individual  horizon  are 
also  reasonably  constant,  field  measurements  of 
horizon  thickness  and  total  mass  of  the  horizon 
per  unit  area  (fuel  loading)  along  with  detailed 
laboratory  measurements  of  particle  density,  sur- 
face area  to  volume  ratio,  and  specific  heat  fully 
describe  the  state  relationships. 

Dynamic  Relationships 

Dynamic  relationships  between  the  state  rela- 
tionships and  flux  properties  are  based  on  analy- 
sis of  the  influence  of  a  particular  state  relation- 
ship on  mass  or  heat  transport  in  soils  when  ex- 
ternal controls  are  held  constant  and  the  flux 
achieves  a  steady  state  equilibrium.  These  are 
theoretical,  idealized  relationships  between  the 
state  relationships  and  the  flux  properties  of  hy- 
draulic conductivity  and  the  void  diffusivities. 
Mass  flux  is  based  on  assumptions  involving  the 
state  relationships  of  porosity,  internal  surface 
area,  and  a  tortuosity  factor  which  describes  ag- 
gregate porosity  microvariations.  Hydraulic  con- 
ductivity of  individual  horizons  is  readily  meas- 
ured. To  relate  these  measurements  to  the  state 
relationships,  the  Kozney  permeability  (Fair  and 
Hatch  1933;  Childs  and  Collis-George  1950) 
model  was  used.  The  intrinsic  permeability  is 


where  "a"  is  an  arbitrary  constant.  Intrinsic  per- 
meability, with  units  of  area,  is  converted  to  hy- 
draulic conductivity,  with  units  of  velocity 
through  the  Darcy  equation  relationship 

K'  =        i^L  (2) 

where  -q  is  the  viscosity,  and  dpldh  is  a  unit  drop 
in  hydraulic  head.  For  water,  K'  =  1  cm/sec  cor- 
responds to  1.02  x  10"6  cm2. 


Similitude  Relationships 

Characterization  of  nonsteady-state  heat  and 
mass  flux  requires  a  time  variable.  The  tradi- 
tional time  variable  used  in  similitude  relation- 
ships is  the  timelag  or  time  response  of  a  system 
to  an  instantaneous  change  from  one  environmen- 
tal state  to  another.  Such  responses  generally  fol- 
low negative  exponential  function  and  are  charac- 
terized by  negative  integer  exponents  of  e  (i.e., 
the  time  to  reach  e"1  of  response  is  the  first  time- 
lag;  the  time  between  the  value  at  e"1  and  e"2  is  the 
second  timelag,  etc.).  Only  rarely  are  the  individ- 
ual timelags  equal  between  the  first  and  last  (e"n) 
because  of  hystereisis  or  differential  responses  to 
external  forces. 

Similitude  relationships  describe  time-depend- 
ent diffusion  processes  in  forest  soils.  The  first, 
the  Fourier  number  (Luikov  1966)  for  heat  trans- 
fer and  the  corresponding  relationship  for  mois- 
ture flux,  is 

Fh=  ~w~  88  °-14  (3a) 

for  heat  transfer,  and 

Fm=  «  0.18  (3b) 

for  moisture  flux,  where  x  is  the  thermal  diffu- 
sivity,  rh  is  the  heat  flux  timelag,  H  is  the  thick- 
ness of  the  layer,  v  is  the  moisture  diffusivity,  and 
7m  is  the  moisture  timelag.  The  similitude  num- 
bers in  equations  3a  and  3b  were  obtained  from 
Fosberg  et  al.  (1970)  and  Fosberg  (1973),  respec- 
tively. These  similitude  relations  describe  bulk 
characteristics  of  the  diffusion  process  but  do  not 
describe  details  of  the  physical  processes  (i.e., 
mass  and  energy  transferred  in  the  voids  as  op- 
posed to  that  transferred  in  the  particles).  These 
relationships  are  general,  however,  since  they  are 
valid  for  bulk  properties  in  both  sorbing  and  non- 
sorbing  porous  media. 

A  second  set  of  similitude  relationships  which 
are  specific  to  sorbing  porous  media,  and  which 
are  more  specifically  related  to  the  physical 
processes  have  been  defined  (Fosberg  1975). 
Structural  properties,  namely  the  porosity,  are  in- 
cluded in  these  similitude  numbers.  Also,  the 
thermal  and  moisture  timelags  of  the  particles 
are  included  in  those  numbers.  These  numbers  re- 
flect the  storage  and  slow  uptake  or  release  of 
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heat  and  moisture  to  the  voids.  The  expression 
for  moisture 

T  V113 

rm      H4/3r^   (1-0)  loc' 

breaks  down  for  nonsorbing  particles  or  when  the 
particle  timelag  becomes  very  short  since  the 
horizon  timelag  is  limited  at  the  low  end  by  mo- 
lecular diffusion  through  the  voids.  Heat  transfer 
for  sorbing  materials  is  similarly 

T  X2'3 

Fh  =  H-Tpi"  (1-*)  <3d) 

where  rm  and  are  the  horizon  timelags  for  mois- 
ture and  heat,  H  is  the  horizon  thickness,  and  rpm 
and  rph  are  the  particle  timelags  for  moisture  and 
heat.  This  relationship  also  breaks  down  when 
heat  conduction  between  particles  becomes  im- 
portant and  the  storage,  uptake,  and  release  of 
heat  becomes  small.  These  last  two  relationships 
most  appropriately  describe  vapor  and  heat 
transport  through  litter,  duff,  and  humus  layers. 
The  Fourier  number  relationships,  while  describ- 
ing all  horizons,  are  the  best  estimators  of  diffu- 
sion processes  in  nonsorbing  horizons. 

To  provide  a  data  base  for  modeling  heat  and 
mass  flow  through  porous  media,  data  from  the 
literature  and  this  study  have  been  combined  to 
calculate  the  flux  properties.  The  basis  for  com- 
bining the  observations  to  obtain  flux  character- 
istics is  defined  by  the  state,  dynamic,  and  simili- 
tude relationships. 


Franco)  on  the  Manitou  Experimental  Forest 
near  Woodland  Park,  Colo,  to  establish  the 
Kozney  constant  in  the  permeability  equation. 
Sample  holders  were  cut  from  standard  snow 
tubes  with  an  inside  diameter  of  5.804  cm  and  a 
length  of  5.230  cm.  One  end  of  each  sample  tube 
was  beveled  to  a  sharp  edge.  By  carefully  trim- 
ming the  duff  and  humus  against  this  edge,  the 
sample  tube  could  drop  through  the  horizon  with- 
out disrupting  the  natural  structure.  Since  the  in- 
tent of  this  experiment  was  to  determine  heat  and 
mass  flux  characteristics,  small  twigs,  bark  plate- 
lets and  other  nonneedle  material  were  not  re- 
moved from  the  sample. 

Tube  depth  (5.230  cm)  somewhat  restricted 
sampling  location  to  those  sites  having  individual 
horizons  of  at  least  this  thickness  since  the  tubes 
must  be  filled  to  reduce  error  in  the  permeability 
measurements.  These  restrictions  biased  the  re- 
sults to  horizons  in  high  basal  area  or  dense 
stands  and  to  near  tree-base  locations.  They  did 
insure,  however,  that  sample  tubes  contained 
only  material  from  a  single  horizon,  and  that  the 
results  could  be  considered  to  be  "homogeneous" 
from  the  standpoint  of  experimental  replication. 

Permeability  measurements  were  made  by 
measuring  the  pressure  drop  between  the  top  and 
bottom  of  the  sample  tube.  Samples  were  placed 
in  a  permeameter  (fig.  1)  which  had  a  mesh  screen 
support  fine  enough  to  prevent  duff  and  humus 
particles  from  being  drawn  into  the  flow  system. 
No  pressure  drop  was  induced  by  the  screen.  Air 
was  drawn  through  the  sample  with  a  vacuum 
pump  at  a  known  flow  rate.  The  pressure  drop 


Experimental  Procedures 

Two  separate  sets  of  experiments  were  carried 
out  to  calibrate  the  dynamic  and  similitude  equa- 
tions previously  described.  The  first  experiment 
was  designed  to  calibrate  the  permeability  equa- 
tion to  make  infiltration  calculations.  The  second 
experiment  was  to  determine  the  thermal  and 
moisture  timelags  of  individual  duff  layers,  and 
was  to  estimate  the  similitude  numbers. 


Permeability  and  Bulk  Density  Measurements 


Exhaust 


Intake 


Sample  holder 


Samples  of  duff  and  humus  layers  were  taken 
from  under  ponderosa  pine  (Pinus  ponderosa 
Laws),  lodgepole  pine  {Pinus  contorta  Doug),  and 
Douglas-fir    (Pseudotsuga    menziesii  (Mirb) 


Figure  1  .—Schematic  ot  permeameter.  Air  is  drawn  through 
sample.  Pressure  drop  is  measured  between  free  air  and 
pressure  at  the  bottom  of  the  sample. 
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Eight  samples  of  Douglas-fir  duff,  seven 
samples  of  Douglas-fir  humus,  eight  samples  of 
ponderosa  pine  duff,  five  samples  of  ponderosa 
pine  humus  and  two  samples  each  of  lodgepole 
duff  and  humus  were  collected. 


Heat  and  Vapor  Flux  Measurements  in  an 
Environmental  Chamber 

To  determine  heat  and  mass  diffusion  rates  in 
duff,  and  to  test  the  hypothesis  that  particle  sorp- 
tion properties  control  the  diffusion  rate  rather 
than  simple  void-controlled  diffusion,  a  two-com- 
partment environmental  chamber  was  con- 
structed and  instrumented  for  temperature,  heat 
flux,  and  humidity.  The  chamber  consisted  of  two 
plenums,  each  approximately  24  cm  deep  with  a 
51  by  41  cm  interface  (fig.  2).  A  ponderosa  pine 


Sensor 
connection 


Figure  2.— Schematic  of  humidity  cabinet.  Sample  of  duff  layer  is  placed  on  screen  in  the  center  of  the  cabinet. 


across  the  sample  was  measured  with  an  EQUI- 
BAR  type  120  pressure  meter.  The  permeability 
measuring  system  was  constructed  so  that  volu- 
metric flow  rate  and  pressure  drop  were  measured 
at  the  downstream  end  of  the  flow  system.  There- 
fore, pressure  fluctuations  induced  by  the  pump 
were  minimal.  Pressure  drops  through  this  open 
flow  were  typically  a  few  tens  of  miUimeters  from 
free  air  pressure.  A  number  of  flow  rates  and  cor- 
responding pressure  drops  were  measured  for 
each  sample.  All  permeability  measurements 
were  made  in  the  morning  to  ehminate  environ- 
mental pressure  fluctuations  from  frequent  after- 
noon thunderstorms.  Intrinsic  permeability  was 
calculated  by  eq.  (2). 

Once  permeability  was  determined  for  each 
sample,  the  container  and  intact  samples  were 
dried  for  24  hr  at  105  °C  to  determine  bulk  den- 
sity. 
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duff  layer  8.7  cm  thick  was  placed  on  a  screen 
shelf  between  the  two  plenums  so  that  diffusion 
would  take  place  only  through  the  duff  sample. 
This  sample  was  lifted  from  the  field  as  an  intact 
horizon.  A  tight  fitting  metal  lip  around  the 
screen  edge  prevented  leakage.  Thermistors  were 
placed  2  cm  above  the  duff  layer  in  the  upper 
plenum;  on  the  upper  surface;  at  distances  2.9, 
3.0,  4.1,  4.9,  and  5.6  cm  from  the  top  surface;  on 
the  bottom  surface;  and  2.2  cm  below  the  sample 
in  the  lower  plenum.  Moisture  content  was  meas- 
ured with  Delmhorst  humidity  probes.  These 
wooden  humidity  sensors  are  0.237  cm  thick  and 
have  a  response  time  of  about  1  hr.  The  moisture 
probes  were  placed  1.9  cm  above  the  top  surface, 
on  the  top  surface,  and  at  distances  2.9,  3.7,  4.3, 
4.6,  and  5.2  cm  from  the  top.  Moisture  probes 
also  were  placed  at  the  lower  surface  of  the  duff 
and  2  cm  below  the  sample  in  the  plenum.  A  heat 
flux  plate  was  placed  3.8  cm  from  the  top  surface. 
All  sensors  were  located  near  the  center  of  the 
sample  to  avoid  any  wall  effects.  The  entire  cabi- 
net was  insulated  with  10-cm- thick  foam  rubber. 

Two  experiments  were  carried  out.  The  first 
was  set  up  to  determine  water  vapor  flux  proper- 
ties of  the  duff  layer,  specifically  the  moisture 
timelag  and  effective  diffusivity  of  the  layer.  Ef- 
fective diffusivity  is  used  here  to  indicate  the 
composite  effect  of  sorption  as  well  as  transport 
processes.  The  duff  sample  was  brought  to  an 
equilibrium  moisture  condition  with  a  saturated 
salt  solution  of  NaC103  at  20  °C  to  produce  a  uni- 
form relative  humidity  of  75%.  Once  the  sample 
attained  moisture  equilibrium,  it  was  held  at  that 
condition  for  48  hr  before  beginning  the  desorp- 
tion  test.  A  saturated  salt  solution  of  CaCl2#6H20 
(equilibrium  relative  humidity  of  32%  at  20 °C) 
was  placed  in  the  lower  plenum,  and  the  NaC103 
salt  solutions  were  removed.  Humidity  was  moni- 
tored on  all  sensors  at  20-min  intervals  for  96  hr. 

The  second  experiment  was  to  determine  the 
thermal  timelag  and  heat  flux  properties  of  the 
duff  sample.  The  sample  was  conditioned  with 
the  NaC103  salt  solution  for  1  week.  After  this 
initial  conditioning,  a  60-watt  light  bulb  in  the 
upper  plenum  was  used  as  a  heat  source.  Tem- 
perature was  monitored  at  2-min  intervals,  and 
heat  flux  was  monitored  continuously  for  6  hr. 

Results 

Intrinsic  permeability  was  calculated  from  the 
pressure  drop  and  flow  rates  with  eq.  (2).  After 
appropriate  unit  conversion 


K=  =  6.936  xlO"7v/Ap  (4) 

where  v  is  in  cm/sec,  Ap  is  in  mm  of  Hg,  and  K  is 
in  cm2.  The  relationship  between  permeability 
and  readily  measured  properties  such  as  bulk  den- 
sity or  other  parameters  used  in  fire  behavior 
models  was  evaluated  through  the  Kozney  formu- 
lations of  permeability.  Measured  bulk  densities 
were  combined  with  data  (Anderson  1969,  Brown 
1970)  on  surface  area  to  volume  ratio  and  density 
of  particles  to  calculate  porosity.  These  trans- 
formed data  were  used  to  calculate  the  calibration 
coefficient  in  Kozney 's  equation  (table  1).  Bulk 
density  of  duff  layers  ranged  from  0.06  to  0.12 
gm/cm3  while  the  permeability  ranged  over  an 
order  of  magnitude  (from  10"3  cm2  to  10"4  cm2).  As 
bulk  density  increased,  permeability  decreased. 
Humus  samples  had  permeabilities  much  lower 
than  duff,  ranging  from  10"5  cm2  to  10"6  cm2.  Also, 
bulk  densities  for  the  humus  sample  ranged  from 
0.1  gm/cm3  to  0.36  gm/cm3  and  showed  a  much 
stronger  density  influence  on  permeability  than 
did  the  duff.  Combining  all  duff  samples  into  a 
single  regression  equation  of  permeability,  versus 
the  group  of  variables  (03/(l-</>)2  a2),  the  slope 
(Kozney  coefficient)  was  2.81  x  10"2  ±  1.5  x  10"2, 
where  the  range  is  expressed  as  a  standard  error 
of  estimate.  The  corresponding  slope  for  humus 
was  4.16  xl0"2±  1.9  xlO"2. 

Diffusion  and  sorption  experiments  in  the  envi- 
ronmental chamber  were  evaluated  from  integral 
moisture  contents  measured  by  the  Delmhorst 
moisture  probes  and  internal  temperature  meas- 
ured with  thermistors.  Both  temperature  and 
moisture  content  were  transformed  to  nondimen- 
sional  variables  of  actual  value  less  the  final  value 
divided  by  the  initial  value  less  the  final  value, 

T-Tf 
To  -  Tf 

and 

M-Mf 
M0-  Mf 

These  ratios,  as  a  function  of  time  (fig.  3  and  4), 
approximated  a  straight  line  in  a  log-linear  graph, 
indicating  that  a  timelag  could  be  determined  for 
both  moisture  and  temperature.  Moisture  timelag 
was  determined  to  average  50  hr  for  the  8.7  cm 
duff  layer  (fig.  3)  and  the  thermal  timelag  was 
found  to  be  1.45  hr.  Calculation  with  eq.  (3a)  and 
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Table  1. — Characteristics  of  duff  and  humus  samples. 


Sampl e 

Bulk 
dens  i  ty 

Pa  rt  i  cl  e 
dens  i  ty 

Surface  area 
and  volume 
rat  io 

Pack i ng      _        .          „         ,  .  , . 
..         Porosity  Permeability 
ratio                    '  ' 
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Constant 

P 

P 

a 

e          <t>  k 

a 

g/orn 

Douglas-fir  duff: 
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cm 


-1 


Dimen- 
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D%men- 
sionless 
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0.6 

69. 1 

0. 163 

0.834 

1  .2 

x 

10  ; 

2.  62 

x 

.098 

.6 

69.1 

.163 

.834 

9-6 

X 

!°: 

2.10 
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.160 

.840 

8.4 
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1.73 
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.  b 

69. 1 

1  J.  0 

.  1  42 
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.  o5o 

1  f 
1 . 6 
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2. 44 
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.118 

.6 

69. 1 

•  197 

.803 
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X 

1 

6.08 

X 

.097 

.6 

69.1 

.162 

.838 

8.2 

X 

;£? 

1.75 

X 

.082 

* 

.6 

69. 1 

•  137 

.863 

1  .2 

X 

IO-4 

1.67 

X 

.  U/3 

.  b 

by.  1 

loo 
.  1 11 

Q-7Q 

.  O  /o 

1  .  b 

X 

i  0 

1     r  L 

1  . 56 

X 
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•  o3 

•  4/ 

1  c\r\ 
1  UU 

.  4yb 

Q  1 

0 .  1 

X 

in"7 

1  cC. 

1 . 56 

X 

or/: 

.  zi>o 
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1  UU 

.  ^55 

1  1 
1  .  1 

X 

o"6 

3-  4/ 

X 

o  eft 

ii7 

1  no 
1  UU 

/1  c  1 

1  c 
1 .  b 

X 

,0-<; 

c  oZ 

5.  26 
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.301 

.47 

100 

.640 

.360 

4.1 

X 

3.60 

X 

.272 

.47 

100 

•  579 

.421 

4.4 

X 

io~6 

1 .98 

X 

humus : 

.  1  /  D 

•  Pb 

1  UU 

1 1  L 

£ft£ 
.  bob 

£  ft 
b .  0 

X 

io"6 

9  nft 
z .  UO 

X 

.  jDD 

1  nn 

1  UU 

.054 

C  9 

j  •  z 

X 

1  u 

C    '5  7 

X 

mus : 

1 1 1 
.mi 

.  b 

JUU 

1  ftc 

ftl  Q 

0.  n 

y .  u 

X 

10"6 

i>  •  by 

X 

1  9  1 

.  b 

.5UU 

.  ZUZ 

7qO 

.  /yo 

7  £ 

/  •  b 

X 

A  1  n 
b .  1  u 

X 

•  119 

.6 

300 

.198 

.802 

2.1 

X 

1 

1 .60 

X 

.121 

.6 

300 

.202 

•  798 

6.3 

X 

5.06 

X 

.  1 06 

.6 
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.177 

.  823 

4. 1 

X 

10 

2.  30 
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.  b 
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X 
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7 
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X 

1  9ft 

.  D 
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X 
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57.6 
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1  .  1 
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X 
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•  51 
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•  796 
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X 
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X 
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.51 

57.6 

•  153 
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X 

10.  J 
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X 
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•  51 
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•  157 
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X 
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X 

.065 
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.127 

.873 
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X 

10  1 

4.26 

X 

.059 

.51 
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.116 

.884 

5.1 

X 
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3.30 

X 

.089 

.51 

57.6 

.175 

.825 

5-1 

X 

!S-i 

9.23 

X 

.077 

•51 
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.151 

.849 
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X 
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X 
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.079 

.56 
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.141 
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X 

-^ 

4.73 

X 
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X 

10 

5.60 

X 
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10 
10 
10 


-2 
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-2 
-2 


10 
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10 


-2 
-2 
-2 
-1 
-1 


10 


-2 
-2 


10 


-2 
-2 

!-* 

-2 


10 
10 
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10 
10 
10 
10 
10 
10 
10 


-2 
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-2 
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2 
3 
3 


10 
10 


-2 
-2 


(3b)  yielded  effective  diffusivities  of  7.5  x  10"5 
cm2/sec  and  2.03  cm2/sec  respectively  for  water 
vapor  and  heat. 

Thermal  conductivity  of  the  duff  layer  was 
measured  near  the  end  of  the  heat  flux  experi- 
ment after  a  reasonably  constant  temperature 


gradient  had  been  established  (a  variation  of  less 
than  0.1°C/cm  hr).  The  thermal  conductivity  was 
found  to  be  2.67  x  IO-5  cal/(cm  sec  C). 

Similitude  relations  expressed  in  eq.  (3c)  and 
(3d)  can  be  used  to  calculate  horizon  timelags. 
Using  horizon  timelags  measured  in  this  experi- 
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Figure  3.— Desorption  of  ponderosa  pine  duff  layer.  Aver- 
age timelag  is  50  hr. 
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Figure  4.  — Normalized  temperature  change  of  ponderosa 
pine  duff.  Average  timelag  is  1.46  hr. 
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ment,  along  with  particle  timelags  for  ponderosa 
duff  particles,  the  similitude  number  for  moisture 
was  1572,  and  for  heat  was  1438.  It  should  be  em- 
phasized that  these  numbers  are  from  a  single  ex- 
periment; therefore,  they  may  not  be  universal. 

Summary 

Hydraulic  conductivities  of  duff  and  humus  are 
large,  ranging  from  a  few  tenths  of  a  centimeter 
per  second  to  as  much  as  30  cm  per  second.  An  at- 
tempt to  predict  these  hydraulic  conductivities 
from  the  readily  measured  property  of  bulk  den- 
sity met  with  marginal  success  in  that  the  stand- 
ard error  of  the  prediction  equation  yields  a  50% 
uncertainty  in  the  hydraulic  conductivity.  While 
the  error  is  large,  so  is  the  conductivity  in  relation 
to  precipitation  rates,  and  the  two  horizons  will 
not  impede  infiltration.  Any  puddling  of  water  in 
these  horizons  must  come  from  water  standing 
above  less  permeable  substrata. 

A  single  pair  of  experiments  on  heat  and  water 
vapor  diffusion  through  a  ponderosa  pine  duff 
layer  support  the  concept  that  heat  and  water 
vapor  transport  rates  through  litter,  duff,  and 
humus  are  regulated  by  the  sorption  properties  of 
the  organic  materials.  The  timelags  which  would 
be  expected  from  these  layers,  if  the  material 
were  perfectly  inert  to  the  heat  and  vapor  content 
of  the  voids,  would  be  very  short.  Reduction  of 
the  diffusivities  according  to  eq.  (2c)  and  (2b)  in 
the  appendix  is  slight.  Observations  of  large  hy- 
draulic conductivities  in  these  horizons  would 
lead  to  the  idea  of  short  timelag.  Measured  time- 
lags  are  much  larger,  and  calculated  effective  dif- 
fusivities are  much  lower  by  several  orders  of 
magnitude  than  would  be  expected  from  an  inert 
solid.  Direct  measurements  of  thermal  conduc- 
tivities also  show  this  reduction  of  heat  flux.  The 
implication  is  that  the  transport  of  heat  and 
water  vapor  by  diffusion  is  controlled  by  the  stor- 
age capabilities  of  the  particles,  and  that  source- 
sink  terms  must  be  included  in  the  traditional  in- 
put-through flow-output  models  for  these  hori- 
zons. 

Results  from  this  limited  heat  and  water  vapor 
diffusion  experiment  qualitatively  support  the 
hypothetical  model.  Little  quantitative  can  be 
said  because  of  the  limited  data.  More  confidence 
can  be  placed  on  the  hydraulic  conductivities, 
however,  because  of  the  substantial  number  of 
samples. 
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APPENDIX 


Basic  Definitions  and  Elementary  Manipulations 

1.  State  Variables— the  properties  are: 

Surface  area  to  volume  ratio  (a) 
Surface  area  of  particle  (S) 
Volume  of  particle  (Vp) 

Particle  density  (qp) 
_  Mass  of  particle  (Mp) 
Volume  of  particle  (Vp) 

Bulk  density  of  horizon  {qb) 
Mass  of  particles  (Mp) 
Volume  of  horizon  (Vb) 

Porosity  of  horizon  (0) 

Volume  of  void  space  (Vv) 
Volume  of  horizon  (Vb) 

Thickness  of  horizon  (H)  =  thickness  of  a  homo- 
geneous layer  in  which  the  above  variables  are 
constant. 


Packing  ratio  ((3) 
Volume  of  particles  in  the  horizon  (EVP) 
Volume  of  horizon  (Vb) 

Multiplying  the  basic  definition  of  porosity  by 
unity  Mp/Mp, 


0  = 


Qp  ~  Qb 


(1) 


2.  Dynamic  Variables 


Void  diffusivities  in  porous  media  (Millington 
1959;  Millington  and  Schearer  1971;  and  Penman 
1940)  are  reduced  from  free  air  values  by  the  ratio 


—  =  —  =  02x 


X0 


(2a) 


These  basic  state  variables  are  interrelated  by  the 
following  definitions: 

Volume  of  void  space  (Vv) 
_       Volume  of  horizon  (Vb) 
~  -  Sum  of  particle  volumes  (EVP) 


where  v  and  x  are  the  vapor  and  thermal  diffusivi- 
ties in  the  voids;  v0  and  x0  are  the  corresponding 
values  in  the  free  air,  0  is  the  state  variables 
porosity,  and  x  is  a  nondimensional  tortuosity 
factor  defined  by  the  implicit  function 

02x-l  +  (l-0)x  =  O  (2b) 


10 


Agriculture— CSU,  Fort  Collins 


-gS  6 
£03  o 

o  n° 
h  ft  • 

ft  «  2 

cd     .  o 

h  *  -3 
go  CO 

es  .-CO 

£  & 

73 

8? 

$  D  «« 


ft 

X 

w 

co 
bo 
3 
03 


«  73 
73  O) 

3  *H 

cd  3 

co  eg 
_  co 

•S  a 

3  55 
73 

3  50 
CO 

mh  bo 
73  P 

O  73 

00  3 
Cp  Co 

23  P 


CO  4h 

73  ^ 

O  cu 

a  J 

.2  «s 

23  co 
ed  o 

«3 

.S  o 

*o  * 

73  hO 

s  5 

CO  O 
CO  M 


00  -u 
CO 

r3  co 

co  o 
co  <4h 

*  M 

CO  CU 

8  % 
£  a 
8 

|£ 

O  73 

c  .a 
s  i 

O  C3* 

CP 


CO 


CO 


£0 


co 

ed 
co 


co 
3 

6 

3 
J3 

73 
C 

ca 


be  0) 
u  'B 

fa 


73 
3 

CO 

u 
O 
fa 


J* 
co 
o 

03 


CO 

3 

a> 

73 

3  co 

73  73 

s  p 

cd 

33  ft 


CO 
co 


cd 

CO 

O 

E| 

^  O 
ft 


O 
ft 

cd 
> 


CO  CO 
53  03 

<4H  73 

CO  S 
CO  Cd 
CO  ^ 

_.  CO 

73  O 
cd  p3 

JS-2 

73  co 

ftQ 


«*H  ft 

CO 

es  c  co 

co  .2  13 
*S  23  O 

2-S  S 
S  g 

„  (S  fl 

O   _,  CO 

■o  eg 

-4J  73 

73  i«  C 
co  73  cd 

'  >»  S 
cd  fl 

co  -2 

Sco  ea 
cj  Si 

CO 


00 

cd 
co 


co 
cd 


u 

CO 

73  S 
CO    >  ^ 

>>  P  -a 


piS  ft  oo 


hi  O 


CO 
O 

CO 

cd 

•  i— i 
o 

CO 

C 

O 

•  i— i 
-u 
cd 
u 


O 

CO 

co" 


CO 

•  p-H 

o 

a 

"co 

00 

73 

o 

CO 


.3  ^> 
Z  & 

$  o 

»-(  ft  - 

a£  a 

tj  .a 

ft  co  O 
03  P3 


P  £ 

oo  73 

73  CO 

cd  3 

+J  CO 

co  cd 

_  CO 

•5  a 

fl  CO 


S3 

cd  . 

*  I 

u  co 

co  CO 


o 


co  53 

73  ^ 

P  CO 

a  J 


.2 

23  oo 

cd  O 

i3  co 

.a  § 


00  -U 
00 

73  co 

oo  o 

CO  «tH 

CO  CO 

CO  ^ 
CO 


J3 
EH 


ft 

s 

p 
o 


cd 


cd  -00 

73  fe  X 

S  CO 

Q  be 

2d  « 
P3 


co 
cd 

u 


3 

J3 

m— i  bfi73  bp 
3  O  (d  p 
co^ 

O  73  3 

60  S  "co 
cp  cd 

23  -  P 
co   C2  B 


O 
ft 

cd 
> 


CO 

s 

3 
J3 

73 

a 

cd 


CO  -fl 
o  8 


73 
S3 
cd 

P 


o 
o 

P3 


a)  u  a)  a> 

3   ft  O 

S    CO  tin  73 

73  I?   09  cd 

3  £     «  +3 

*a      73  o 


^  cd  ^ 
^  oo  cd 
co  o 


S3  g  ®  .  53 
CO  'P    >  tC 

ft  ftQ 


l^1 

O  73 
'  CO 
.3 
S3  3 

o  a1 

CD 

co  2; 

«4H  ft 

co 

cd  c  co 

_»  .2  "g 

2  js  a 

_Q  'CO 

Si 

-U    Cd  . 

73  J?,  3 
CO  73  cd 

'  >^  ft 

5-i  co 
cd  3 
£  * 

CO  ^ 

H  co  cd 
CO  ,o 

ed  S3  co 
^  co  +i 

55  73  « 
CO    >  ^ 

o  ™ 


CO 

3 

CO 


3 

CO 

cd 
co 


co 
co 
3 

JS 
Is 

X! 

U 

CO 
-u 
cd 

JJ 

•i— i 
P 

CO 

3' 
O 
-u 
cd 


o 

00 

co" 
u 
3 

CO 


CO 

3 

CO 

73 
tw 
O 

CO 


ft  . 

*-<  S3 
ft  co  CJ 


3 


O 


73  CO 

3  a 

cd  3 

-t-»  GO 

oo  cd 

.a  s 

3  ^ 


oo  53 
73  73 

O  CO 

a  I 

3  a 
O 


cd 
u 


CO 

o 

.  CO 

3  5 
ft 


oo 

-u>  co 

"3  co 

3  s- 

co  o 

CO  <+H 

CO  CO 

2  ft 


J3 
EH 


a 

o 

CO 


> 

^   CO  • 

co  CO  3 
ed  -00 

c3 

73 


3  _  ft  3 

*r  ^3  « 


ft 

S  co 
g  CO  3 

K  M  P3 


3  «  o 

«w  bD73  bp 
"3333 
©  cd  O 


l^1 

O  73 

'  CO 

3 
3 
U* 

co 

3  .& 


O  73 
3 


cd  » 

►  o 

co  n 
3 


t> 


1  § 

co  "3 
-°  o 


oo  73 

a  * 

3 


ft'oS  co  •§  ^  So 


o 
fa 


CO 

o 
03 

ft 
o 


o 
a 

cd 

M   CJ   S3  > 

CO   ft  CO  CO 

^   CO  P,  "3 
r-H  **h  cd 

31  ^    »H  ^ 

-3  ft  O 

a  co  <3  73 

if  co  3 

73  1?    K  Cd 

3  J2  «3  +j 

>»  g  cd  2 
33  "a  cd  ch 

«  co  «  3 


fa  S  p  J3 
ft  ftQ 


73  2 

w  J2 
cd  3  co 

co  .2  V, 
£  43  P 

cd  cd  3 

a  co 

3  ^ 

cd 

O    C  « 

^  cd  3 
O  co 

+3    03  _. 

co  w  cd 
3  ^  S 

S  3  3 

3  co  cd 

CO  ,Q 

co  T; 

3  co 
co  -3 

73  « 

co  > 
cd  P 

i— H      ft  CO 


CO 

cd 


CO 
CO 

3 

JS 

Id 
u 

CO 

cd 

j* 

•i—i 
P 

CO 

3" 
O 

•fH 

cd 


P 

00 

co" 

Sh 

3 

-u 

oo 


CO 

3 
fa 

00* 

73 
u 
O 


CO 


.a  »^ 

^  c^> 

CO  ■ 

ft 
o  r 

u  ft 
ft  cd 
4J  fa  .3 

in  S3 

cd 


DO 


E: 

u  co 
co  CO 


cd 

73 

3 

cd 


u 
o 
fa 

< 
P 

%™ 


CO 


M 
fa 
CO 


CO 


05  3 

3 


bo 
3 
cd 
03 

73 
3 

cd 

Sh 

O 
fa 


^  P 

«  73 

73  CO 

3  ^ 

cd  3 

-U>  CO 

co  cd 

_  CO 

.a  a 

a  % 

3  ^ 


73 

S3  3 
cd  J2 
vm  bo 

P  73 

co  a 

Cd 


GO  <H 

73  ^ 

O  cu 

a  I 

3  a 

.2  «3 

23  GO 

cd  O 

23  co 

a  § 


ro-3 
73  bo 
3  3 
cd  P 

CO  j3 
co  O 

a  ft 


cd 


CO 

cd 

bfi.co 


Ja3 
o 
o 

03 


33  CO 

2  3 
co  ft 
^  a) 

a  co 

-°  be 

73  73 

3  © 

cd 

>>  s 

^3  ft 


-  O 


cd 


GO 

■+J  CO 

73  co 
3  s-i 

GO  O 
CO  <+h 

^  X 

CO  CO 

£  a 

8 

P  73 

'  CO 

,3 
3 

cr 

CO 

VM  ft 

CO 

GO  ^ 

cd  f3  oo 

ft  .s*g 


S3 

o 

co 


^2 

cd 
co 


co  '3  ft 
2  o  t— i 


cd 

CO 

a  2 

co  ^3 
fa  3 

ft 


a  > 

%  CO 

3 

«C  cd 
&  T3 

co  a 

go  cd 
cd  ^ 

&  Id 

co 

-2-3 
O 

cd  c 

73  co 

>  «a 

O  «*H 

D.3 


_Q  CO 

E3 

Id 

o  3  co 

_  cd  3 

o  -  co 

+j    Cd  . 

73  i«  3 

CD  73  cd 

H  >»  ft 

S-H  CO 

j  cd  3 

co  cd 

S    CO  Jh 

3  co  cd 

CO  ^ 

CO     2  >H 

cd  S3  co 

£  CO 

^  73  g 

^3 
l^ftg 


3 

co 
cd 


co 

CO 

3 

JS 

"cd 

CO 

cd 

•  I-H 
P 

GO 

3" 
O 

•I-H 

cd 


P 

00 

co" 

3 

-u 

GO 


CO 

3 
fa 

co 
73 

Si 

O 


CO 


